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Abstract 

A Monte Carlo algorithm for calculating the single-particle spin-echo 
small-angle neutron scattering (SESANS) correlation function is presented. 
It is argued that the algorithm provides a general and efficient way of cal- 
culating SESANS data for any given shape and structure. 
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1 Introduction 



Spin-echo small-angle neutron scattering (SESANS) has recently emerged as a 
new way of applying neutron scattering to the investigation of the structure of 
matter [?, ?]. The method is particularly useful for large structures in the size 
range from 10 nm up to several microns. This is the same size range covered 
by techniques like light scattering and ultra-small anglc^ nciitron scattering (US- 
ANS). The use of neutron spin echo in measuring elastic scattering, however, 
renders beam coUimation unnecessary, thus avoiding the low fluxes from which 
USANS suffers. In comparison to light scattering, the use of neutron allows 
for study of opaque or highly concentrated samples. The SESANS method is 
presently on an active developing stage, and the theoretical concepts and meth- 
ods from which conventional small-angle neutron scattering (SANS) benefits 
have only recently started to be derived [?] and be applied to the analysis of ex- 
perimental data [?]. As in the case of SANS, the data analysis can be performed 
with model-dependent or model-independent methods. In the latter case, the 
experimentally obtained scattering functions are inverted to obtain a curve rep- 
resenting the pair distance distribution function p{r). The most well-known 
realisation for this procedure is the Indirect Fourier Transform by Glatter [?]. 
Similar analysis can be carried out by the maximum entropy method [?] and 
the regularization method of Svergun [?]. Model- independent analysis is par- 
ticularly simple in the case of SESANS, as the relation between the SESANS 
correlation function G{z) and the small-angle scattering correlation function 
7(r) is given by an Abc^l integral equation [?], for which there exist standard 
numerical methods for solution. 

In model-dependent analysis, mathematical functions that model the scat- 
tering intensity from a system of particles with presumed shape, structure, and 
ordering are fitted to the experimental data. From the fitted parameters one 
obtains information such as the size and shape of the particles, their inner struc- 
ture and size distribution, and the inter-particle interactions that create ordered 
structures. In SANS terminology, the last piece of information is contained in 
the structure factor, whereas information pertaining to single-particle scattering 
is contained in the form factor. Analytical or semi-analytical functions for the 
scattering form factors or scattering amplitudes are known for several geome- 
tries [?] . These functions are easy to extend to include core-shell structures and 
polydisperse assemblies. 

As will be elaborated in the following section, G{z) is related to the SANS 
scattering cross section (d) / (d) by a two-dimensional cosine transform [?] . Know- 
ing the full detectable (d)/(d) as an analytical function or as tabulated values of 
(d)/(d) vs. Q, calculation of G{z) is straightforward by numerical integration. 
In principle, it would be desirable and more efScient with analytical expressions 
for G{z) for different geometries, analogous to form factors for SANS, so that 
SESANS data can be analyzed with a similar tool-box of model functions. An- 
alytical expressions for G{z) for scattering from single homogenous and hollow 
spheres have already been derived using the concept of the mass correlation 
function 7(r) [?]. If 7(r) of a given structure is known, calculation of G{z) is 
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simple. For non-spherical geometries, however, expressions for 7(r) can take 
complicated forms [?, ?, ?]. Moreover, it is not an easy task to extend j{r) 
from homogenous to multi-domain structures. 

In this paper we investigate an alternative method to calculate G{z). It is 
for all practical purposes of general validity and straightforward to implement 
for any shape and structure. The idea is to perform a Monte Carlo calcula- 
tion of the pair distance distribution function p(r), from which the SESANS 
correlation function G{z) and also the SANS scattering cross section (d)/(d) 
can be calculated by a single numerical integration. Monte Carlo methods have 
been used by several authors in calculation of SANS spectra [?, ?, ?, ?]. The 
algorithm for calculating SESANS curves is outlined in Section 3. In Section 4 
we present results of the calculations for different shapes and structures. 



2 SESANS Theory 

The measured quantity in a SESANS experiment is the loss of polarization a 
neutron beam suffers by being scattered by the sample. By passing through 
magnetic fields before and after interacting with the sample, the neutrons are 
subjected to Larmour spin precession. In the case of no interaction with the 
sample, the precessions before and after the sample area cancel each other, 
yielding a spin echo preserving the polarization state of the beam [?]. The 
presence of a scattering sample produces precession lengths differences that are 
functions of the scattering angle. The ensuing depolarization is a function of 
the SESANS correlation function G{z) [?]: 

^ = exp [G(z) - G(0)] = exp [G(0) (Go(z) - 1)] , (1) 

-TO 

where Go{z) is the normalized correlation function. The relation between the 
SANS macroscopic scattering cross section (d)/(d) and the SESANS correlation 
function has already been derived [?, ?]: 

where A and t denote the wavelength of the neutron beam and the thickness 
of the sample, respectively. Qy and Qz are the cartesian components of the 
scattering vector Q, the incident beam lying along the x axis. The integrations 
in ^ are in practice defined by the area in the yz-plane covered by the detector. 
The spin-echo length z is a function of the neutron wavelength, the sample 
position, and the configuration of the magnetic fields [?] . 

We now consider a system of non-interacting particles isotropically embed- 
ded in a homogenous matrix or dispersed a solvent. The SANS scattering cross 
section can be written in terms of an intra-particle form factor P{Q) and an 
inter-particle structure factor S{Q): 

^^(Q)=n,V^P{Q)S'{Q), (3) 
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where is the number density of scattering particles and V is the volume of 
a particle. Most of the analytical structure factors have been calculated for 
systems of monodisperse spheres. The effective structure factor S'{Q) includes 
approximate corrections to S{Q) due to particle polydispcrsity or anisotropy 
[?, ?, ?]. We will in the following consider dilute system, for which we may 
ignore inter-particle scattering and set S'{Q) = 1, corresponding to an ideal 
gas. The form factor P{Q) is related to the average size and shape of the 
individual particles and to their inner structure. Focusing on a single particle, 
the form factor can be written in terms of the density correlation function 7(r): 

P{Q) = l%{rf-^4nr'{d), (4) 

where 7(r) is defined by [?, ?]: 

lir) = \ ^W^P (r') Ap (r' + r)^ , (5) 

where the braces () denote averaging over all orientations of the position vector 
r. Ap(r) is the scattering length density at a position r inside the particle, 
minus the constant scattering length of the surrounding medium (in most cases 
a solvent). D is the largest chord length of the particle, so that 7(r) = for 
r > D. For a homogenous particle j{r) is proportional to the overlap volume 
between the particle and its identical "ghost" that has been shifted by a distance 
r. For an inhomogenous particle, the volume of the overlapping region must be 
weighted with the product of the scattering length densities of the respective 
regions [?]. An important identity is [?]: 

{ApfV= / 7(r)47rr2(d), (6) 
Jo 

where Ap is the difference between the average scattering length density of the 
particle and that of the surrounding medium (in most cases a solvent). The 
normalized density autocorrelation function 70 (r) is defined through 7(r) = 
(Ap)^7o(r) and has the property 7o(0) = 1. 

In this paper we focus on the single-particle contribution to the SESANS 
spectrum. A direct real-space interpretation of G{z) was presented by Krouglov 
et al. Tim2003. For a system of non-interacting particles, the SESANS cor- 
relation function can be written G{z) = X'^tnpGp{z), where the single-particle 
SESANS correlation function bears the following relation to the structure func- 
tion 7(r): 

/oo 
7 Ux^+zA (d) 
-00 

= 2{ApfvJ^ Jo{^/x^+^){d). (7) 



4 



Note that both G{z) and Gp{z) are dimensionless. Knowing j{r), the SANS 
form factor and the single-particle SESANS correlation function can be calcu- 
lated. An important quantity in SESANS is the total scattering probability, 
given as G{0) [?, ?]. It relates to the observed depolarization and thus gives an 
indication on the magnitude and detectability of the SESANS signal. From the 
above we have: 

G(0) = 2X^tnpV / 7(r) (d) 
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= 2X't<t>iApr / 7oM(d) (8) 
Jo 

= Xhcj^iApfl (9) 

where (f) — UpV is the volume fraction of the particles and I is the mean length 
of all chords contained in the particle [?]. Finally in this section, we remark 
that in the case of an ensemble of polydisperse, non-interacting particles, the 
equations above take the form 

/ .VD^ \ 

G(z) = 2XHnpUApfv 70 ( V^^TI^ j (d) \ (10) 

G(0) = 2XH^^(^{ApfvJ%oix)id)y (11) 

where () now stands for the averaging over the particle sizes. The number 
density is given by np ^(l)/{V). lf{Apy is the same for all particles. Equation 
lfTT]l can be written G(0) = A^i0 (Ap) 1^, where Jw is the weight-averaged mean 
chord length. 



3 Monte Carlo calculation of G{z) 

3.1 Calculating the pair distance distribution function 

The pair correlation function 7(7') is related to the pair distance distribution 
fimction (pddf) p{r) by 

p{r)=r^j{r). (12) 

p{r) is the probability of two random points within the particle being separated 
by a distance r. By random sampling of distances between points within the 
particle, and keeping statistics of the sampled distances, p{r) can be found for 
any particle shape. The interval < r < D is partitioned into M + 1 histogram 
bins, indexed from to M, D being the maximum distance between two points 
belonging to the particle. M — 200 was used for the calculations presented in 
this paper. 

Geometrical points are sampled uniformly from a volume that circumscribes 
the volume of the given particle as closely as possible. When Nm = 1000 points 
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have been selected, the points that fall outside the shape function of the particle 
are discarded, leaving Nr points. The points are sampled in batches of in 
order not to exhaust the computer memory. The shape of the sampling volume 
is essential, as it ensures that most of the sampled points will belong to the 
particle. This makes the calculations far more efficient than sampling from a 
circumscribing rectangular box, discarding the points that fall outside the shape 
function of the particle [?, ?]. The sampling volumes are rectangular, spherical, 
or cylindrical, depending on the shape of the particle. Uniform sampling from 
a given distribution or volume by an inverse method is a well-known technique 
in Monte Carlo calculations [?]. For completeness, we give a brief outline of the 
procedure in Appendix^ For the spherical and cylindrical sampling volumes, it 
is straightforward to limit the samling to given segments or sectors. Regardless 
of the shape of the simulation box, it is necessary that it encloses the entire 
particle ensure that different parts of the particle contribute to p{r) in proportion 
to their volume. This can also be accomplished by allocating a fraction Vi/V of 
the randomly selected points to each domain i. This alternative is particularly 
suited for multidomain structures with inhomogenous density distributions. 
The algorithms proceeds by calculating the Nr{Nr — l)/2 distances defined 

by the accepted points. The distance d = \J \Ya — rb|^ between the points 
and Tf, is counted into the histogram, wheighted by the product of the scattering 
length densities pa and pb of the two points. The update of the histogram is 
carried out by the following algorithm: 

p{i)^p(^i^mt(^^+Q.5^^+PaPb, (13) 

where the int() function represents truncation of decimals, leaving the integer 
index i of the histogram bin. The index i is related to the intraparticle distance 
by r = iD/M . Adding 0.5 to the argument of the int() eliminates the problem 
overcounting smaller distances because of the decimal truncation. 

When the Nr{Nr — l)/2 distances have been counted, the program checks 
the accuracy of the calculated p{i) after a criterion to be described below. If 
the accuracy is not accepted, an new batch of Nr points is selected, and the 
resulting new Nr{Nr — 1) distances are counted into the histogram following 
equation 1)13(1 . The total number of sampled points belonging to the particle, 
ie accepted points, is denoted Np; and the total number of sampled points is 
denoted Nm- 

The normalization of p{i) is done so that the form factor will satisfy P{Q = 
0) = 1. This is accomplished by calculating p(i) —p{i)/ (CD/M), where 

C=J2 PaPb = {ApfN„iNp~l)/2 (14) 

a=l b=a+l 

is the sum of the distance weights. The last equality in H14|l is valid in the limit 
of large Np. For a homogenous particle, C will be proportional to the number 
of sampled intra-particle distances. Knowing C, (Ap)^ can thus be calculated. 
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In most cases the volume V of the particle will be known beforehand. If not, it 
can be found from the MC calculations by the relation V — NpVhox/NM, where 
Vbox is the volume of the simulation box. 

Numerical tests showed that a reliable test for the accuracy of p{i) was 
to compare p{i)N, the pair distance distribution function calculated from A'' 
sampled points, with p{i) n-iooq- This was done for every time p{i) had been 
calculated with 1000 new points and upgraded with lOOOx (1000— 1)/2 distances. 
The calculations were halted when 



n:^^(')N-P(')N-loo^^ <0.25%. (15) 
V T,i [pCOat-iooo] 

For homogenous particles, the required number of points lay around Np — 5000. 
For inhomogenous particles, the number could be significantly higher; for core- 
shell particles with equal volumes and opposite signs of the scattering length 
densities, values up to Np — 40000 were typical. Nevertheless, calculation of a 
full data set Go{z) took only a few seconds on a notebook equipped with a 750 
MHz Pentium-Ill processor. 



3.2 Finding 7(r) 

To find j{r) from the calculated p(r), Eq. 1)12(1 faces us with the problem of 
dividing by zero or small values of . We overcome this problem by exploiting 
the small-r expansion of 7(r): 

7(r) = a + br + cr^ -(-... (16) 

The expansion parameters are known as differential parameters and are related 
to structural features of the particle [?]. The second order parameter c will 
in most cases be zero, and the last term could thus be neglected or replaced 
with a third-order term dr^ . For the sake of the generality of the algorithm, 
however, we apply the expansion as given in equation (|16|l . One should bear in 
mind, though, that there are particle shapes for which an expansion like H16|l 
can not be carried out [?]. However, it is valid for most realistic geometries. We 
find the differental parameters a, 6, c by polynomial least-squares fitting [?] of 
ar^ _l_ j,j,3 _|_ ^jjg Monte Carlo calculated p{r) at small r. For the sum of 

squared residuals we have: 

1=0 

where = iD/M . is minimized with respect to a, b, and c; and the resulting 
linear system is solved for a, b, c with standard routines [?]. The summation in 
H17|) runs from i — to i = m, where the index number to < Af is decremented 
until the sum of squared residuals Em(a, b) has a sufficiently low value. When 
this is the case, 7(r) up to index z = to is given by a+bri+crf . For the remaining 
values of r, 7(r) is calculated directly by p(r)/r'^ . An example is demonstrated 



7 



in Figure ^ In the calculations in this paper, the initial value of m is set at 
half the distance between r = and the first peak of p{r). Finally, to facilitate 
a consistent computation of G{z) and G(0), jo{f) = l{r)/a is calculated by 
normalization. 



3.3 Calculation of G'(O) and G'o(^) 

The single-particle SESANS correlation function Gp{z) is calculated by numer- 
ical evaluation of the integral lO, neglecting the prefactor 2{Apf'V. Having 
found 7o(r), we need to evaluate 70 (V + z^) for arbitrary values of x and 
z. This is accomplished by natural cubic spline interpolation [?]. The integral 
lO is evaluated using a 16-point Gauss-Legendre quadrature [?, ?]. The nor- 
malization of Gp{z) is carried out by dividing by Gp(0), which is calculated at 
the beginning by evaluating for z = 0. Subsequent calculations of Gp{z) 
are normahzed by dividing by Gp(0), thus yielding Go{z). The total scattering 
probabihty G(0) is given by 2 (Ap)^ npX'^tGp{0), where the calculation of (Ap)^, 
if unknown a priori, can be done by means of Equation 1141) . 



4 Results and discussion 

In Figure 121 we have plotted the MC-calculated Go{z) curve for a homogenous 
sphere with radius R = 50 together with the analytical expression for Go{z) [?]. 
The latter reads: 



Go(C) = 




where ^ = z/R. The agreement between the calculated data and the analyti- 
cal function is excellent. The other curves represent MC-calculated Go{z) for 
spheres with R = 50, but with radial density profiles following a hyperbolic form 
p(r) — r^". The density profile has a pronounced effect on the appearance of 
the Go{z) curves and also on the total scattering probability G(0). For the full 
sphere, one has I = 3R/2 = 75, whereas for p{r) — r~^, the calculations gave 
1= 5.4. 

Figure|3|shows the MC-calculated Go{z) for hollow spheres with outer radius 
i? = 50 and varying inner radii. Comparing Go(z) of the hollow spheres with 
that of the full sphere, the interesting feature is the appearance of a small 
shoulder at a z-value corresponding to the inner diameter of the hollow spheres. 

Go{z) for core-shell spheres with inner radius Ri = 30 and outer radius 
i?2 = 50 are plotted in Figure 0] for different combinations of the scattering 
length densities pi and p2 ■ Some values of pi , p2 give rise to strong oscillations 
in Go(z), which were also found in the multishell calculations by Rekveld et al. 
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MTR2000. The quantitative interpretation of these oscillations was given by 
Uca [?]. Minima arise from correlations between particle regions with opposite 
signs of their scattering length densities. Correspondingly, maxima are related 
to identical signs in different regions, or to correlations within the same region. 
The positions of these extrema give information on the typical distances between 
or within these regions. Thus, the position of minimum at z = 40 in Figure 0] 
is related to the typical shell-core distance, which is taken as the core radius 
plus half the shell thickness, which in the case shown in the Figure is exactly 
i?i + (i?2 - Ri) /2 = 40. The maximum at 2; « 70 is due to the shell-shell 
correlations, of which the typical distance, taken as the core diameter plus twice 
the half shell thickness, is 80. The core-core correlations are expected to give a 
maximum around z = i?i = 30, but this is hidden by the inital part of Go{z). 

Alternatively, oscillations in Go{z) can be interpreted in terms of the dif- 
ferential parameters introduced in Section 13.21 Inserting the expansion Hlfi|l 
into (0 shows that each term in H16|l . except the zero-order term a, produces a 
contribution to Go (z) containing a maximum. The position such a maximum is 
shifted to higher z for higher order terms. Thus the minimum of Gq{z) at z = 40 
in Figure ^ can be attributed to a large negative differential parameter of high 
order in the expansion Accordingly, the maximum at 2 = 70 comes from 

a positive differential parameter of even higher order. Relating such high-order 
differential parameters directly to structural features of the particle in a unique 
way is, however, a challenging task in small-angle scattering theory [?, ?, ?]. 

MC-calculated Go(z) for homogenous and hollow cyhnders are plotted in Fig- 
ure It should be noted that for cylinders with a high aspect ratio L/2R, the 
maximum of p(r) is shifted towards small r. As a consequence, the least-squares 
method for determining j(r) described in Section ^01 must use a small number 
of data points. However, this is not found to pose a serious problem. The Go{z) 
curve for the homogenous cylinder shows the same characteristic features as in 
the model calculations by Uca et al. Uca2003. For the hollow cylinder, however, 
there are remarkable differences. The SESANS functions are concave at small z, 
and there is a clear shoulder at a position corresponding to the inner diameter. 
For the thinnest of the cylindrical shells, the shoulder appears to give a disconti- 
nous first derivative at z = 100. In the corresponding p{r) curve, the maximum 
also appeared discontinous. Because of this, it was necessary to increase the ac- 
curacy of the numerical quadrature described in Section 13.31 to avoid numerical 
artefacts. At the highest z-values, the curves practically coincide, all showing 
the loss of correlations characteristic for anisotropic structures. 

SESANS functions for core-shell cylinders analogous to the core shell spheres 
in Figure 01 are plotted in Figure The inner and outer radii are i?i — 30 and 
i?2 = 50, and the cylinder length is i = 250. The positions of minima and 
maxima can be interpreted in the same way as for the spheres. For the lowest 
curve, the correlations are almost completely lost already at z = 100. This 
can be related to the fact that for this particular combination of core and shell 
volumes and scattering length densities, the overall scattering length densities 
of the core and shell are nearly equal, but of different signs. Thus for z greater 
than 100, corresponding to correlations only along the cylinder axis, the core 
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and shell contributions cancel each other out. 

In Figure |3 the SESANS functions for various triaxial ellipsoids are shown. 
The case for ellipsoids of revolution have been discussed by Uca et al. Uca2003. 
The SESANS functions of the ellipsoids with small eccentricities have been 
calculated by sampling points from a circumscribing sphere. For higher eccen- 
tricities, a circumscribing cylinder was used. Sampling random points from a 
confocal ellipsoidal coordinate system [?] would in principle eliminate the prob- 
lem of sampling and discarding points that fall outside the ellipsoid. However, 
applying the inversion method of Appendix ^ to such special coordinate sys- 
tems would in most cases require the numerical solution of nonlinear equations, 
thus loosing the advantage of sampling efficiency. 

We close this section with a discussion of the efficiency of the Monte Carlo 
method, with the prospective application of analyzing real SESANS data. The 
calculations in this work were performed using the random number generation 
ran3 given by Press et al. NumericalRecipes. This is a very reliable routine, 
but for the purpose of rapid calculations, simpler generators could be applied. 
Although the Monte Carlo method is general, the examples shown are for rela- 
tively simple structures. Additional refinements are possible, but to the cost of 
increased computing time. Polydispcrsity would require multiple calculations of 
7(r), as given by Equation H1(J|) and Hll|l . but this would in principle be required 
to be done only once for each G'o(z)-spectrum. The presence of radial density 
profiles p{r) could be included already in the Monte Carlo sampling procedure, 
as described in Appendix IXI For a sphere, this method requires that r^p(r) is 
described by an analytically invertible function. If not, an additional sampling 
must be performed, preferrably from a distribution that closely follows p{r), to 
obtain the desired density distribution. Again, this calculation needs to be done 
only once for each p{r). 

5 Conclusion 

Knowledge of analytical expressions for the SESANS correlation function Gq{z) 
for a host of geometries and structures would be ideal and efficient for calculating 
model curves and analyzing experimental SESANS data. As such expressions 
are hard to obtain for complex structures due to the difficulty of deriving gen- 
eral expressions for 7(r) to be used in equation |(7J), one has had to apply the 
scattering functions from conventional small-angle scattering to equation ||2Jl 
obtain Gq{z). The Monte Carlo algorithm outlined in this paper represents an 
alternative method, which is general and straightforward to implement from the 
shape function of the geometry in question. The method does not need any of 
the special functions that frequently comes with the SANS scattering functions 
contained in (d) / (d) . 
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A Uniform sampling by the inversion method 



Standard random number generators provide the user with random real num- 
bers uniformly distributed between and 1 [?]. Wishing to sample from a 
distribution function f{x) defined on or limited to the interval a < x < b, one 
sets /(x)(d) — (d) and integrates, getting the following relation between the 
cumulative distribution F(x) and the random variable ^: 

Ja /(a;')(d)' 

Solving the inverse equation x = F~-^{£^), one can convert the random variable 
< ^ < 1 to random variables x uniformly distributed in f{x). 

Considering a sphere with radius R, the three spherical coordinates r, (f>, 9 
have the probability distributions 3r^/R^, l/27r, isin6', respectively. A three- 
dimensional uniformly distributed random variable ^ is then converted points 
uniformly distributed throughout the volume of the sphere by 

= 2^6 (20) 

cos 6* = 1 - 2^3 

When the geometry in question is, say, a hemisphere, a spherical sector or a 
spherical shell defined by the coordinates (i?i,$i,9i) and (i?2,^2,62), the 
distribution functions can be limited to these regions. Equation (|20|l then gen- 
eralizes to 

</> = 6($2-*i) + $i (21) 
cosO = cos 01 — ^3 (cos 01 — cos 82) . 

A further generalization is the inclusion of a radial distribution of scattering 
length by a function p(r). Taking the function p(r) = r~", where a < 3, 
which in particular applies to spherical star polymers and polymeric micelles 
[?], the radial coordinate r follows the probability distribution 47rr^p(r)/F = 
(3 - a) r2-"/i?3-a^ yielding in the case of a hollow sphere 

r-= [^i(i?3-"-i?3-a)]i/(3-")_ (22) 

In the case of a cylinder with radius R and length L, the coordinates r, </>, z 
follow the respective distributions 2r/R^, l/27r, l/L, giving the relations: 

4> = 2^2 (23) 

z = 6^ 

The corresponding generalizations to a cylindrical subvolume are: 

r ^ [i,{Rl-Rl)+Rl]"' 

= 6($2-*i) + $i (24) 
z = 6(^2-^1) + ^!. 
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For a rectangular box defined by a lower corner (ai, a2, as) and an upper corner 
(61,62,63), the transform relations for the cartesian coordinates {xi,X2,X3) takes 
the simple form 

Xi = {bi - tti) + ai. (25) 
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Figure 1: Right: Monte Carlo calculated pair distance distribution function 
p{r) for a homogenous sphere with radius R = 50. The solid line represents the 
polynomial [a + br + cr^) that has been fitted to p{r) up to r = 26.5. Left: 
density correlation function 7(r) calculated from the p{r) curve. The low-r part 
is calculated from the polynomial fit. 
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Figure 2: Monte Carlo calculated SESANS correlation function Go{z) for a 
homogenous sphere with radius i? = 50 (circles). The solid line represents the 

analytical expression for Ga{z). The other curves represent spherical particles 
with a radial density distribution p{r) = r~" and a maximum radius R = 50. 
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Figure 3: Monte Carlo calculated SESANS correlation functions for one ho- 
mogenous (open circles) and various hollow spheres. 
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Figure 4: Monte Carlo calculated SESANS correlation functions for a core-shell 
sphere with inner radius Ri = 30 and outer radius R2 = 50 with different 
scattering length densities. 
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Figure 5: Monte Carlo calculated SESANS correlation functions for one ho- 
mogenous and several hollow cylinders with outer radius R2 = 50, different 
inner radii -Ri, and length L = 250. 
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Figure 7: Monte Carlo calculated SESANS correlation functions for a various 
triaxial ellipsoids. 
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